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We shall revisit the conventional adiabatic or Markov approximation, showing its intrinsic failure 
in describing the proper quantum-mechanical evolution of a generic subsystem interacting with its 
environment. In particular, we shall show that — contrary to the semiclassical case — the Markov 
limit does not preserve the positive-definite character of the corresponding density matrix, thus 
leading to highly non-physical results. To overcome this problem, we shall propose an alternative 
adiabatic procedure which (i) in the semiclassical limit reduces to the standard Fermi's golden 
rule, and (ii) describes a genuine Limblad evolution, thus providing a reliable/robust treatment of 
energy-dissipation and dephasing processes. 

PACS numbers: 03.65.Yz, 72.10.Bg, 85.35.-p 



Present-day technology pushes device dimensions to- 
ward limits where the traditional semiclassical or Boltz- 
mann theory [l[ can no longer be applied, and more 
rigorous quantum-kinetic approaches are imperative Q. 
However, in spite of the quantum-mechanical nature 
of electron and photon dynamics in the core region of 
typical solid-state nanodevices — e.g., superlattices Q 
and quantum-dot structures 0| — the overall behavior 
of such quantum systems is often governed by a com- 
plex interplay between phase coherence and energy re- 
laxation/dephasing [Bj|, the latter being also due to the 
presence of spatial boundaries Therefore, a proper 
treatment of such novel nanoscale devices requires a the- 
oretical modeling able to properly account for both co- 
herent and incoherent — i.e., phase-breaking — processes 
on the same footing. 

The wide family of so-called solid-state quantum 
devices can be schematically divided into two main 
classes: (i) a first one which comprises low-dimensional 
nanostructures whose electro-optical response may be 
safely treated within the semiclassical picture 0| (e.g., 
quantum-cascade lasers Q), and (ii) a second one 
grouping solid-state devices characterized by a genuine 
quantum-mechanical behavior of their electronic subsys- 
tem (e.g., solid-state quantum logic gates [9|) whose 
quantum evolution is only weakly disturbed by decoher- 
ence processes. 

For purely atomic and/or photonic quantum logic 
gates, decoherence phenomena are successfully described 
via adiabatic-decoupling procedures in terms of 

extremely simplified models via phenomenological pa- 
rameters; within such effective treatments, the main 
goal/requirement is to identify a suitable form of the 
Liouville superoperator, able to ensure/maintain the 
positive-definite character of the corresponding density- 
matrix operator [Tl| . This is usually accomplished by 
identifying proper Lindblad-like decoherence superoper- 
ators [ll|, LL2J , expressed in terms of a few crucial system- 
environment coupling parameters [l3j ]. 

In contrast, solid-state devices arc often characterized 



by a complex many-electron quantum evolution, result- 
ing in a non-trivial interplay between coherent dynam- 
ics and energy-relaxation/decoherence processes; it fol- 
lows that for a quantitative description of such coher- 
ence/dissipation coupling the latter need to be treated 
via fully microscopic models. 

To this aim, motivated by thepower and flexibility of 
the semiclassical kinetic theory [l[ in describing a large 
variety of interaction mechanisms, a quantum general- 
ization of the standard Boltzmann collision operator has 
been proposed Q; the latter, obtained via the conven- 
tional Markov limit, describes the evolution of the re- 
duced density matrix in terms of in- and out-scattering 
superoperators. However, contrary to the semiclassical 
case, such collision superoperator does not preserve the 
positive-definite character of the density-matrix opera- 
tor. 

To overcome this serious limitation, in this Letter we 
shall propose an alternative adiabatic procedure which (i) 
in the semiclassical limit reduces to the standard Fermi's 
golden rule [lj|, and (ii) describes a genuine Limblad 
evolution, thus providing a reliable/robust treatment of 
energy-dissipation and dephasing processes. 

In order to discuss the main features and intrinsic limi- 
tations of the conventional adiabatic or Markov limit, let 
us recall its general derivation following the fully operato- 
rial approach proposed in Given a generic physical 
quantity A — described by the operator A — its quantum 

plus statistical average value is given by A = tr|j4p|, 
where p is the so-called density-matrix operator. Its time 
evolution is dictated by the total (system plus environ- 
ment) Hamiltonian. Within the usual interaction picture, 
the latter can be regarded as the sum of a noninteracting 
(system plus environment) contribution plus a system- 
environment coupling term: H = H a + H'\ the corre- 
sponding equation of motion for the density-matrix oper- 
ator — also known as Liouville-von Neumann equation — 
in the interaction picture is given by: 



dp 1 
~dt 
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2 



where TL denotes the interaction Hamiltonian H' written 
in units of h. 

The key idea beyond any perturbation approach is that 
the effect of the interaction Hamiltonian & is "small" 
compared to the free evolution dictated by the noninter- 
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acting Hamiltonian H a . Following this spirit, by formally 
integrating Eq. ([1]) from — oo to the current time t, and 
inserting such formal solution for p l (t) on the right-hand 
side of Eq. |T]), we obtain an integro-differcntial equation 
of the form: 



dt' 



W{t), H*(t'),?(t') 



(2) 
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We stress that so far no approximation has been intro- 
duced: Equations fl} and ([2|) are fully equivalent, we 
have just isolated the first-order contribution from the 
full time evolution in Eq. ([I]). 

Let us now focus on the time integral in Eq. ©. Here, 
the two quantities to be integrated over t' are the in- 
teraction Hamiltonian TL 1 and the density-matrix op- 
erator p 1 . In the spirit of the perturbation approach 
previously recalled, the time variation of p l can be 
considered adiabatically slow compared to that of the 
Hamiltonian TL written in the interaction picture, i.e., 
Wit') = U 1 (t')nU (t')- indeed, the latter exhibits rapid 
oscillations due to the noninteracting evolution operator 

U Q (t) — e b 2- . As a result, the density-matrix operator 
p l can be taken out of the time integral and evaluated 
at the current time t. Following such prescription, the 
second-order contribution to the system dynamics writ- 
ten in the Schrodinger picture for the case of a time- 
independent interaction Hamiltonian TL comes out to be: 



dp 
~dt 



TL, 



)C,p 



(3) 



with 
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dt'TL l (t') = 2 / dt'U {t')TLUl{t') . (4) 



As discussed extensively in [16| , the operator AC describes 
energy-conserving scattering events (real processes) as 
well as energy-renormalization contributions (virtual pro- 
cesses); the latter are known to play a minor role, and 
in general may be safely neglected; such approximation 
amounts to impose the following time-reversal symmetry 
on the system dynamics: TL l (t) = TL l (—t) 17]. Within 
such approximation scheme, the operator /C in |4| may 
be rewritten extending the time integration from — oo to 
+oo: 



K 



+ 00 



dt'U {t')TLUt{t') 



(5) 



The effective equation in ([3]) has still the double- 
commutator structure in ^ but it is now local in time. 



The Markov limit recalled so far leads to significant 
modifications to the system dynamics: while the exact 
quantum-mechanical evolution in (TTJ) corresponds to a 
fully reversible and isoentropic unitary transformation, 
the instantaneous double-commutator structure in © 
describes, in general, a non-reversible (i.e., non unitary) 
dynamics characterized by energy dissipation and de- 
phasing. However, since any effective Liouville super- 
operator should describe correctly the time evolution of 
p and since the latter, by definition, needs to be trace- 
invariant and positive-definite at any time, it is imper- 
ative to determine if the Markov superoperator in Q 
fulfills this two basic requirements. As far as the first is- 
sue is concerned, in view of its commutator structure, it 
is easy to show that this effective superoperator is indeed 
trace-preserving. In contrast, as discussed extensively in 
[IH , the latter does not ensure that for any initial condi- 
tion the density-matrix operator will be positive-definite 
at any time. This is by far the most severe limitation of 
the conventional Markov approximation. 

By denoting with {|A)} the eigenstates of the noninter- 
acting Hamiltonian H , the effective equation ([3]) written 
in this basis is of the form: 



dp Xi\ 2 _ 1 
dt ~ 2 



[PxiX 2 ,x' 1 x 2 Px , 1 x 2 - 'Px 1 x 2 ,x' 1 x> 1 px 2 x 2 ]+H-v- 

X[X' 2 

(6) 



with generalized scattering rates given by: 

'Px 1 x 2 ,x' 1 x' 2 = -rH' XiX ^H'^ 2X , 2 8(£\. 2 - e X '. 2 



(7) 



e\ denoting the energy corresponding to the noninteract- 
ing state I A). 

The well-known semiclassical or Boltzmann theory [l[ 
can be easily derived from the quantum-transport for- 
mulation presented so far, by introducing the so-called 
diagonal or semiclassical approximation. The latter cor- 
responds to neglecting all non-diagonal density-matrix 
elements (and therefore any quantum-mechanical phase 
coherence between the generic states Ai and A2), i.e., 
PAiA 2 = /ai^AiA 2I where the diagonal elements fx de- 
scribe the semiclassical distribution function over our 
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noninteracting basis states. Within such approximation 
scheme, the quantum-transport equation ([6]) reduces to 
the well-known Boltzmann equation: 



Markov superoperator in Q cannot be employed, since 
it does not preserve the positive-definite character of the 
density matrix p\ x \ 2 . 



dfx 
dt 



= X)( p aa'/a'-Pa'a/a) 



where 



Pxw = V 



2tt 

aa.a'A' = -r\H'\\\ <5 (ca - £A' 



(8) 



(9) 



are the conventional semiclassical scattering rates given 
by the well-known Fermi's golden rule [IBj . 

At this point it is crucial to stress that, contrary to 
the non-diagonal density-matrix description previously 
introduced, the Markov limit combined with the semi- 
classical or diagonal approximation ensures that at any 
time t our semiclassical distribution function f\ is al- 
ways positive-definite. This explains the "robustness" 
of the Boltzmann transport equation ([8]), and its ex- 
tensive application in solid-state-device modeling as well 
as in many other areas, where quantum effects play a 
very minor role. In contrast, in order to investigate gen- 
uine quantum-mechanical phenomena, the conventional 



As anticipated, aim of the present Letter is to propose 
an alternative formulation of the standard Markov limit, 
able to provide a Lindblad-like scattering superoperator, 
thus preserving the positive-definite character of our den- 
sity matrix. To this end, let us go back to the integro- 
differential equation ((2]). As previously discussed, the 
crucial step in the standard derivation is to replace /S l (i') 
with p l (t). Indeed, since in the adiabatic limit the time 
variation of the density matrix within the interaction pic- 
ture is negligible, the latter can be evaluated at any time. 
Based on this remark, what we propose is the follow- 
ing time symmetrization: given the two times t' and t, 
we shall introduce the "average" or "macroscopic" time 
T — ^1*1 and the "relative" time r = t - t' . The basic 
idea is that the relevant time characterizing/describing 
our effective system evolution is the macroscopic time T. 
Following this spirit, it is easy to rewrite the second-order 
contribution in Eq. Q in terms of the new time variables 
T and r: 



J 



dr 



T -7 



(10) 



r 



Following again the spirit of the adiabatic decoupling, we 
shall now replace p 1 (T — At) with p l (T); the resulting 
effective equation rewritten in the original Schrodinger 
picture comes out to be: 



dT 



p{T) 



dr 



(ii) 

Neg lecting again renormalization contributions (see note 
[17l|). Eq. (fTl") may be rewritten by extending the time 
integration over r from — oo to +oo: 



dr 



W I -r 



-^r),p(T) 



(12) 

By Fourier expanding the above symmetric convolution 
integral we finally get the desired Lindblad-like scattering 
superoperator: 

1 



with 

t{u>) 



dp 
dT 



/27T 



dio 



drU n 



C(u>), C{bj),p 



(13) 



(14) 



We stress how the proposed time symmetrization gives 
rise to a fully symmetric Lindblad-like superoperator (ex- 
pressed in terms of the operator C only), compared to the 
strongly asymmetric Markov superoperator in ([3]). 

If we now rewrite the new Markov superoperator in 
(|12p in our noninteracting basis A, we obtain again the 
effective equation of motion in ([6]) , but now the general- 
ized scattering rates in (0 are replaced by the following 
symmetrized quantum scattering rates: 



£Ai + £ a 2 e K + e K 



(15) 

The above scattering superoperator can be regarded as 
the quantum-mechanical generalization of the conven- 
tional Fermi's golden rule; indeed, in the semiclassical 
limit (Ai = A2,A' 1 = X' 2 ) the standard formula in is 
readily recovered. 

At this point a few comments are in order. As dis- 
cussed extensively in 16[ , also for the simplest case of a 
standard two- level system — i.e., a generic quantum bit — 
the standard Markov superoperator predicts a non-trivial 
coupling between level population and polarization de- 
scribed by the so-called T3 contributions. In contrast, for 
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a two- level system coupled to its environment, the pro- 
posed quantum Fermi's golden rule does not predict any 
T 3 coupling term, thus providing a rigorous derivation of 
the well-known and successfully employed TiT 2 dephas- 
ing model Moreover, it is imperative to stress that 
in the presence of a strong system-environment interac- 
tion the adiabatic decoupling investigated so far needs to 
be replaced by more realistic treatments, expressed via 
non-Markovian integro-differential equations of motion 
(i.e., with "memory effects") Again, while for purely 
atomic and/or photonic systems it is possible to iden- 
tify effective non-Markovian evolution operators 14J , for 
solid-state quantum devices this is still an open problem. 

To summarize, we have critically reviewed the stan- 
dard adiabatic or Markov procedure, showing its intrin- 
sic failure in describing the proper quantum-mechanical 
evolution of a generic subsystem interacting with its en- 
vironment. More specifically, we have shown that within 
the Markov approximation the density-matrix operator 
is not necessarily positive-definite, thus leading to highly 
non-physical results. To overcome this serious limita- 
tion, we have identified an alternative adiabatic proce- 
dure which (i) in the semiclassical limit reduces to the 
standard Fermi's golden rule, and (ii) describes a gen- 
uine Limblad evolution, thus providing a reliable/robust 
treatment of energy-dissipation and dephasing in state- 
of-the-art quantum devices. 

We are grateful to David Taj for stimulating and fruit- 
ful discussions. 
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